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A new class of static plane symmetric solution of Einstein field equation generated by a perfect 
fluid source is put forward. A special family of this new solution is investigated in detail. The 
constraints on the parameters by different energy conditions are studied. The classical stability of 
this solution is discussed. The junction conditions matching to Minkowski metric and Taub metric 
are analyzed respectively. 
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I. INTRODUCTION 

Q ■ Though considerable amounts of vacuum solutions of Einstein equation are known, but, frankly speaking, the 
physical interpretation of many of them remains unsettled pj . As emphasized in :2] , the key to physical interpretation 
^ is to find out the nature of the sources which generate these vacuum spaces. Singularity theorem is one of the 
most important and profound theorem in classical general relativity. For a deeper understanding and control of the 
singularities, the physical characters, such as mass densities, internal pressures, of the sources are imperative. For 
example in black hole solutions, where only singularity "generates" the whole space, investigating of how a matter 
distribution gives rise to the black hole space is necessary to judge the physical significance (or lack of it) of the 
complete analytic extensions of these solutions. 

Concerning how we obtain a comprehending of the sources, there is really no alternative for constructing an interior 
' solution for a matter distribution which matches to the vacuum space in question. We will present a new class of 
04 I plane symmetric solution sourced by perfect fluid. One may be curious about where we can find a realistic plane 
1 symmetric space. In the researches of the early universe, people find there may be several-time phase transitions when 
' the universe is cooled down, and many topological defects, such as monopole, cosmic string, and domain wall may 
OO , survive today. Theoretically, equivalence principle, as principle of general relativity, says that a gravity field can be 
' simulated by an accelerating effect. For the case of a uniformly accelerating observer, an exact simulation should be 
plane symmetric (but curved) space. 

The static vacuum plane symmetric solution, Taub solution was obtained more than 50 years ago [3| . The solution 
is asymptotically fiat but there is a serious time-like singularity in Taub space. People believe that the singularity 
should be superseded by some matter source, for some attempts to seek the sources, see 0- Taub solution permits 
a three-dimensional Euclidean symmetry group, but it is not homogeneous. Here a solution is called homogeneous if 
it permits three translational space-like Killing vectors. Taub solution only permits two: It is not homogeneous in z 
direction. Further, it is shown that generally speaking any singularity free source with refiective symmetry for plane 
symmetric vacuum space must violate dominant energy condition 

In electromagnetics the electronic field of an infinite homogeneously charged fiat board is homogeneous. By analogy, 
in linearized gravity theory the solution for a (5-function like source has been derived 5] , which is spatially homogeneous 
and seems to contradict with 4] . It is suggested that the counterpart of the Einstein theory to the gravitational field 
of a massive Newtonian plane should be described by the Rindler solution rather than Taub's 6]. 

We will present a new class of plane symmetric solution with perfect fiuid source, which definitely violates the 
dominant energy condition. We study some various properties of this new solution and do some preliminary studies 
on how to glue it to a vacuum space. 
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This paper is organized as follow: In the next section we will present our solution. Because the energy momentum 
tensor is very miry and hence the physical properties of this solution are difficult to analyse, the full form of the 
energy momentum tensor is casted into the appendix. We then give a special class of the solution as a starting point 
of the further discussions. In section III we study the constraint on the parameters by the different energy conditions. 
In section IV, we investigate the stability property of this solution. In section V, we try to glue this sourced solution 
to vacuum solutions. The cases of Minkowski and Taub are discussed respectively. A summary and some discussions 
are presented in section VI. 



II. THE SOLUTION 



We find that the following plane symmetric metric 



az^ ^ y-S-L-dh{z) 



{dx^ + dy^) , (1) 
solves the Einstein field equation with a perfect fluid source. Here a, w, c, d are 4 real parameters, and 

/(z) = ce-'^" + ^«e'^^ (2) 

/i(z) = -ce^"^ + u;e°^ (3) 

Clearly there are four Killing fields '§x^ Sij ^"^^ ~V'§x^^^- '^^^ latter three span a Euclidean group G3, which 
implies the plane symmetry. The metric ([TJ solves the Einstein field equation with a source in perfect fluid form 

T^{p{z)+p{z))U®U+p{z)g, (4) 

where T denotes the energy momentum tensor of the fluid, U stands for four-velocity of the fluid and g denotes the 
metric tensor. Since the exact forms of p{z) and p{z) are very involved, we just present them in the appendix. The 
energy momentum tensor of metric ([l} is too complicated to make further analysis. Alternatively we may concentrate 
on a special family of solution ([1]). 

First, we set it; = 0, and hence ([1]) reduces to 

where we have used 



h{z) arctan (e°^ 

lim --^ -A ^JLA = (6) 



Then, we rescale the coordinates t, x and y by 

ct t, ecx X, ecy y. (7) 

So ^ further reduces to 

ds^ = -e-^-^dt" + + g-2(a.+de— ) (^^2 _j_ ^y2y (g) 

Introducing new variables 

a' = -a, y = -d, (9) 

and for convenience labeling a' by a and b' by b, we derive a special solution as follow, 

ds^ = ^e^-^dt^ + dz^ + e^'-'^'+'^^'Hdx^ + dy^). (10) 

There are two free parameters in the above solution. With the above metric p{z) and p{z), respectively given by (|74p . 
(I75p for the original metric ((!]) , are reduced to extraordinarily simple form 

p{z) = -0^(3 + 8fee"" + 36262"^), (11) 
piz) = a2(3 + 46e''^ + 62e2'^^), (12) 

where we have set SttG = 1, and G, as usual, denotes the Newton gravitational constant. In this article we will 
concentrate on metric ([TU|) . Before studying the detailed properties of pH]) . we first investigate the various limitations 
of it. Obviously when a = pI7|) degenerates to Minkowski metric. When = 0, it becomes 

ds^ ^dz^ + e^''''{-dt^ +dx^ +dy^), (13) 

which is just anti-de Sitter (AdS) metric in an unattractive coordinate system. Surely it describes a homogeneous 
and isotropic solution in all directions, not only a plane symmetric one. 
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III. ENERGY CONDITIONS 



More and more cosmological evidences imply that exotic matters violating energy conditions, such as weak, strong, 
dominant or null energy conditions, do exist in our universe. And the partition of the exotic matter takes a much more 
value than that of dust matter. But different requirements on the energy-momentum tensor have been used to prove 
several pivotal theorems, for example singularity theorem and positive energy theorem, in classical general relativity. 
It may be helpful to investigate the properties of the matter in the source by different energy conditions. Further, it is 
pointed out that the source of Taub must violate the dominant energy condition if the source is reflective symmetric 
and without spacetime singularity [3| . Now we begin to study the parameter regions permitted by weak, strong and 
dominant energy conditions. Note that weak and strong energy conditions are two independent conditions. In this 
section we consider a plane source with finite thickness, which resides in the region z > 0. Our solution is an interior 
solution of plane source, thus the pressure p{z) should vanish naturally at some distance from the "ground", z = zq, 
that is, we should match it to a vacuum solution. This condition implies 



p{zo) = a^(3 + 46e"^° + b^e^"'") = 0. 
The two roots of the above equation reads 



bi 



-3e- 



(14) 



(15) 
(16) 



bi and 62 denote the two branches of the solution which have finite thickness. We call them "little branch" and "large 
branch" , respectively. 

To show the physical meanings of the parameters, we calculate the mass per area of this slab a. 



a 



where p{z) is defined in pT|) . For the little branch b = ~e the integral yields. 



For the large branch, the result is 

« = Ire""'" 

72 



-5e 



^2+3az„ _^ g-2e (l4e°^0 2e^"^" + 6^°^" - 12) 



191e 



(17) 



(18) 



(19) 



We see that the mass per area is proportional to a, so it can be treated as a mass parameter of the slab. Whether 
the stuff filled in the slab is ordinary or exotic in some sense depends on the value of azg. Different energy conditions 
impose different constraints on azQ. 

Weak energy condition (WEC) requires that the energy momentum tensor of the source T satisfies, 



T{Z, Z)>0, 



for any time-like vector Z, which equals 



p = -a (3 + »oe 
p + P= -a^i^be"' + 2b^e^''') > 



(20) 



(21) 
(22) 



First, we consider the little branch — 
we investigate three cases depending on the sign of a. 
Case I: a > 0. Inequality ipTj) requires 



We then replace 6 by — e in (PT|) and (P^ . For further discussions 



< 



4-77' 

and inequality (|22p is satisfied spontaneously in this settlement. 
Case II: a < 0. Inequality ((2T|) requires 



(23) 



4+V7' 



(24) 
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Inequality ([22)) requires 

e"^" > 1/2, (25) 

therefore the permitted range of a and zq is 

ga^o > 1/2. (26) 

The case a = degenerates to Minkowski space, which is a trivial case satisfying any energy conditions. In the 
interval, 

e-o e [1/2, (27) 

WEC always can be satisfied for any real a. 

Second, we consider the large branch b = — 3e~"^". Mimicking the above discussions for the case of little branch we 
find that WEC can not be satisfied for any a and zq. 

Then we turn to the strong energy condition (SEC), which requires 

Ric{Z,Z)>0, (28) 

where Ric denotes the Ricci tensor of metric (llOp . and Z is an arbitrary time-like vector. The condition (j28p equals 
to 

/9 + p= -a2(46e"^ + 262e2«^) > 0, (29) 

p + 3p = a^{6 + Abe''^)>0. (30) 

First, similarly, we consider the little branch b = —e^°-^°. Also, we investigate the three cases depending on the 
sign of a, respectively. 

Case I: a > 0. The inequality has been discussed in the case WEC, which is satisfied naturally. One finds ([501) 
is also satisfied naturally. So SEC imposes no constraint on the parameters in this case. 
Case II: a < 0. Inequality ([29]) has been analyzed before, which requires 



^azo > 1/2, (31) 

and inequality ([30[) requires 

gazo > 2/3. (32) 

The case a = degenerates to Minkowski space, as we have pointed out, which is a trivial case marginally satisfying 
any energy conditions. For arbitrary a the requirement of SEC always can be satisfied in the interval, 

e°^° G [2/3,oo). (33) 

Second, we consider the branch b = —3e~°'^° for SEC. Mimicking the above discussions we find that SEC can not be 
satisfied for any a and zq. We see the case b = —3e~'^^° may violate both WEC and SEC for any values of parameters 
a and zq. 

The dominant energy condition (DEC) requires 

P>\pI (34) 

For removing the calculation of absolute value, we consider the cases of p < and p > 0, respectively. First, p < 
requires that 

(& + e-"^'')(& + 3e-°^°) < 0, (35) 

whose solution is an empty set. Second, for the case p > 0, DEC requires 

3 + ebe"' + 2b^e^''' < 0, (36) 

which does not have non-empty solution set either. Hence we conclude that this solution always violates DEC. 
Recalling that [4] demonstrates that any reasonable source of Taub solution must violate DEC, one sees that our 
present solution just satisfies this requirement. 
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IV. STABILITY 



In this section we give an initiative discussion on the global stability of the solution. The basic idea of our 
investigation is best presented by considering an analogy. A Newtonian celestial body, for instance, the Globe, is 
stable, which can be understood as follow: If the Globe undergoes an small perturbation, expanding or contracting a 
little, the total energy (including gravitational energy and the energy of the particles of which the Globe consists) of 
new configuration is higher than the original configuration. Hence the Globe is stable in view of least energy principle. 

A similar argument can be applied to our solution. The first problem is the energy of the gravitational field. 
People recognized that there is no covariant energy momentum of gravitational field long ago. Different expressions 
correspond to different boundary conditions in the Hamiltonian approach of quasi-local mass [sj . Here we adopt the 
Landau-Lifschitz expression Q. According to Landau-Lifschitz expression, the energy momentum of gravitational 
field reads 



1 r 



/npa yp _ pa pp _ pa pp U„M/3 7 _ „Mf „/37^ i „M/3 „7a ('p'' yP J_ F"^ VP ~ —V'^ VP \ 

l^J^/37-Lap J^/3p^7a '-f3a'-ip)\9 9 9 9 )+9 9 {'^ Pp'^ -^a + -fa'^ f3p ap^ /S-y ^ apj 



iP /37 



1//3 7Q/pP pp I pp p- 



P 

7""^ (ip 



pp pP 



^ pf^ ap) + 9 9 y^pa^ip ^ p^^ ap) 



(37) 



where T denotes the affine connection of metric (jlOp and 5^1^ denotes the component of (jlOp . The non-zero components 
of r are 



pt 



zt ~ tt 



2az pz _ _ pz 



-agx 



rL = rL = a(i-h6e-), 

Substituting ([38]) into ([37|) we derive the energy density of gravity, 

Pgra--5° = a2(4 + 46e'^^) 
Then, in the little branch the density per area of gravity is 



r^, = a(l + 6e-). 



Performing the calculation directly, we obtain [llj 

P=2a (e"2+2azo 



which is also proportional to a. For the large branch, the area density becomes 



which yields. 

The total mass associated with the unit slab is a — a + (3. In the little branch, by using and (HI]), 



(38) 
(39) 

(40) 
(41) 
(42) 
(43) 



a = a + f3 



lie 



— 2+3azo 



12 



(44) 



Imitating the case of a Newtonian star, we impose a perturbation to zqj 
^2 



5a 



-2-2e"''^0-2a2o 



^_-|^2g2 ^ ^p2+a2o C2^^^1aza _j_ 2g2+3a2o _|_ g2+4azo _j_ He^"^ °^o+4azo^ j 

Numerical calculations show that stagnation points dwell at zq = — 0.292a^^ and zq = — 1.245a^^. Further we find 

<0, 



(45) 



l5^0'|ro = -0.292a- 
'5^0'|ro = -1.245a-i > 0, 



(46) 
(47) 
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where we have used o < 0, since zq > 0. We conclude that the configuration described by ([10]) is stable when the 
parameters a, tq are constrained by arg — —1.245. The resulting stable metric is a one parameter family, 

In the large branch, similar to the little branch, the total mass of a unit slab reads. 



a = a + P 



(49) 



72 

Imposing a perturbation to zq, we reach 

2 

6a = " g-6-6e-"^o_2azo fQj2e^ + 108e^+°^° + 18e^+2''^« + 6e^+3'''° + e^+^''^° + 3356^"'""+^''^°) Szq. (50) 

Numerical calculations show that two stagnation points for Sa = inhabit at zq = — 0.3163a"^ and zq = 0.9770a~^. 
Further, one can show 

5^cr|r„ = _0.3163a-i > 0, (51) 
3^'^\ro=0.9770a-^ > 0, (52) 

therefore, both of the two stagnation points are stable. They correspond to two spaces of one parameter family, 

rf4ai = -e^'^'dt^ + dz^ + e2(— 4-ii6e"^)(dx2 + dy'), (53) 

respectively. 

V. MATCHING TO VACUUM SOLUTIONS 

Surely, ^ and (fTO|) are rigorous solutions with perfect fluid sources, which can be filled in the whole space. However, 
there is a true singularity if interval of z is not confined. In fact, the Ricci scalar R reads, 

i?= -2a2(6 + 106e"^ + 36^e^""). (55) 

When z goes to oo (—00), Ricci scalar will be divergent for a positive (negative) a. We hereby consider the case that 
this solution is only valid in a finite region and the spacetime is vacuum out of this region. 

The gravitational field must satisfy two boundary conditions: 1. The metric is continuous across the boundary 
surface, and 2. The extrinsic curvatures measured by the different sides of the boundary surface relate to each other 

by 



[K - htr{K)]^ = -T, (56) 

dz ^ dz 



in which K denotes the extrinsic curvature of the boundary, h = g — -S- ® represents the induced metric on the 



boundary, r is the energy- momentum tensor confined to the boundary, and [ ]^ denotes the jump at the boundary, 
i.e., for a quantity Q, [Q\^ = lim(2_2o)^o+ Q{^) ~ l™(z-2o)^o- Qi.^)- The exterior vacuum space has to be plane 
symmetric to match to metric (jlOp . First we study the most simple case that the vacuum out of the source region is 
Minkowskian geometry. Then we analyse the junction condition matching to the well known non-flat plane symmetric 
space: (static) Taub space. 

Before discussing the boundary condition between the slab and the vacuum, we impose a Z2 (reflective) boundary 
condition at z — 0. The continuity condition is naturally satisfied and the jump condition generates 

[X-Mr(if)r ^-l(r)°, (57) 

where 0-1- labels the value of a quantity at z = going from the positive direction, and denotes the value of a 
quantity at z = 0. Using the above equation we derive 

r; = 2diag(2a(l + &), a(2 + 6), a(2 + &)), (58) 
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where b = — e or b = — 3e "^f, depending on the little or large branch. Now we write the Minkowskian metric in 
the following chart, 

dslji^ = -inn?dt^ + dz^ + n'^{dx^ + dy'^), (59) 
where m and n are positive constants. The continuous condition yields 



m — e 



(60) 

which corresponds to "time contraction" and "length contraction" , though to, n, may be greater than 1 so that they 
represent "time dilation" and "length dilation" . From the jump condition ([55]) we obtain t in induced chart by ([55]) 
in the little branch, 

r'; = diag(0, -a, -a), (62) 

while in the large branch, 

t""^ = diag(4a, a, a). (63) 

T vanishes when a = 0, which is consistent with our expectation because the interior metric degenerates to Minkowski 
when a = 0. 

Next we discuss the conditions matching to the Taub metric. The Taub's metric reads, 

ds^ ^ -z-^/^k^dt^ + dz^ + z'^l^l^ {dx^ + dy^) , (64) 
where k, I are two positive constants. The continuous condition requires 

3azo 1 + 21n/cZ, (65) 

in the little branch, and 

2 

azQ — 1 + - In fc^, (66) 
o 

in the large branch. 

The jump condition gives, in the induced chart by ([M]) for the little branch 

= diag ^2(7^ + a - ^z^^'^. + « ^ ' (^''^ 

where zq is given by (|65p. And for the large branch 

r-; = (a+izo-i)diag(-4, -1,-1), (68) 

where zq takes the value in (|66p . If we require the matching is perfect, that is, boundary energy momentum vanishes, 
we arrive at azo — —1/3. 

On the other hand, if perfect matching condition is imposed in the little branch, one needs an infinitely thick slab, 
which can be derived from ([S7|). 

VI. CONCLUSIONS AND DISCUSSIONS 

In this article we present a global plane symmetric solution of Einstein field equation with a perfect fluid source, 
which we interpret as the source of some plane symmetric vacuum space. Four Killing vectors, including a time like 
Killing vector, are permitted in this solution. We find a chart in which the metric is written in time coordinate 
orthogonal form. Then we give a special and much simpler class of the original solution as the starting point for 
further investigations. 

We find the ranges of the parameters in which WEC and SEC can be satisfied, respectively. Interestingly, we find 
that DEC is always violated, no matter what values the parameters are taken. 
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We study the stability properties of this solution. For the little branch, the stable point dwells at zq = —1.245a~^, 
and for the large branch, the stable points inhabit at Zq = — 0.3163a~^ and Zq = 0.9770a^^. 

We do some primary researches on matching to vacuum solutions. Minkowski's and Taub's are studied respectively. 
Furthermore, for clarifying the physical contents of this solution, we should study the weak field approximation to 
find some thing like Newton force among test particles. 

The general form of the plane symmetric solution llj with a perfect fluid source must imply more rich physical 
and mathematical properties. It may deserve to study further in the future. It should be noted that in the present 
article we have examined only the classical stability under z-direction perturbations, though it is the most important 
perturbations of plane symmetric space. The other forms of classical perturbations, quantum and semi-classical 
perturbations should be thoroughly investigated further. 
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VII. APPENDIX 



One derives the Einstein tensor of metric whose components read. 



m'^c^du^ + Adb^&'d^u^ 



A5b^c^d^u^ + 18bc% 



,5„,2 



' + 126^/2^5/2^5^ _^ ISb^c^d^- 
b^Vbcdu^^v + m^/^c^l^u''v - 



le&^c^e^- _ l2Q{bcy/^du\ + Gb'cu^ + 6b''/\'/^u\ - 2Ab^c^du ~ Sbc'/u + 2Ab^/^c'>/^v + 
kic^u-^v - Ub'^/^c'^^'^u^v + Abb^c^u^v^ + Sb^c^u^ + 2'lb'^c^du'' + GOb'^c'^d^u'^ - 90b^/'^c'^^^dv + 
b^l'^c'I'^du + ^b'^cdu^ - S&c^ - eVtec^w-^w + \%b^cu^v^ + Zc^vT'^v^ + Gbti' - Wb'^c^u^ + 



cu 



Sbc'^d^u"^ - 1863/2c9/2„^ - 6b'^^c''^^u-\ 

V- 

+8b^c'^u^ - 36b''cVbcdu''v - SbWbcv + 606^ c^u"?;'^ - Qb''cdu^ a''{bc)-'{c + bu^) 



m^c^u^ + 366wi%2 _ 36^3/2^11/2^ _^ m'^c'Vb'cU^v - IQb^c^u^v 

(,3 3„,4„,2 



Ab^c^u^ - m^c^du^ + Qbc^d/u 

5„2„,7 



'dQb'^c^Vb'cdu^v ~ 2Ab^c^u'v + 



(69) 



Abc^du-^ - Qb^c^d^ - 2b^Vb^u'^v + m^c^Vb^dv/u - Sb^c^u^ + 2y/b'cc^ u'"^ v - Qb^e'd/u- 



Abc^ 



Gb'^c^du^ ~ 15b^c^u^v^ 



Sb'^c^du^ 



2b^cdu' + 2b^Vb?Au^v + Ub^cVbcu^v 

„5„ -2„,2 l6„,8„,2 



Gb'^cVbcu^v - Ab^cu^ + Vbce' 



!.3„4 



du - Ibb^c^v^ 



c 1,6 2 72 6 
DO C d U 



ebc^u-'^v^ - b''u''v^ - Ab^cdu^ + SOb^c^Vbcdv + Ab'^c^du^ + Sb^c^u^ + Ab^cu^ + Qbc^Wcv/u + 



Ube'y/bcu-^v - Ibb^e'd^ + iOb'^c^Wcu^v + I 



!,3 2 



^/bcu^v + Abc^\/bcu ^ - 2bc^Vbcdu 



b'^cd^u^ - Ab^c^du + Ab^Vb?M'^v + 12b^c^d/u + Ab^c^VTcuv - bc^ d 



20b*c*d^u^ 



,7^2^-4 



where u = e"^, v — arctan u 



20b^c^u'^v^ 



and 



b^c^ 



bcuv — A&'Vbcu 



l2bc^Vbcv/u + 
db'^cu^v^] a2(6c)-i(c + 5u2)-4 ^ 



~" y ^ z 



(70) 



(71) 



The Einstein tensor is checked by software Maple. We see that, though the Einstein tensor is rather messy, it reveals 
a significant property of metric ([1]): it describes a solution of plane symmetric sourced by perfect fluid, whose energy 
momentum is in the form Then in coordinate system ([1]), 



Gt, 



(72) 
(73) 



rpx _ rpy _ rp 

X y : 

We use the convention SttG = 1 throughout this article. So, in the comoving system of the perfect fluid, U = d/dt, 

P 



-T* 



P 



(74) 
(75) 
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^ ^ -e^-^dt" + e^'^^'+''^"\dx'' + dy^), 

on the two covers, respectively. Note that the induced metric is not well defined on the joint lines gluing the slab and the 
two covers and the jointing lines between the walls. Hereby the boundary energy momentum tensor is not well defined 
there. Just as well the six lines are zero-measure sets, and then we can omit the energy associated with it without changing 
the result. With this construction it is easy to show that the quasi- local gravitational energy enclosed by the slab and is 
exactly the same as the result by using Landau-Lifschitz expression for our settlement. 



